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Abstract 

The 2iV-dimensional quantum problem of N particles (e.g. electrons) with 
interaction /3/r 2 in a two-dimensional parabolic potential ujq (e.g. quantum 
dot) and magnetic field B, reduces exactly to solving a (2N — 4)-dimensional 
problem which is independent of B and uq. An exact, infinite set of relative 
mode excitations are obtained for any N. The N = 3 problem reduces to that 
of a ficticious particle in a two-dimensional, non-linear potential of strength j3, 
subject to a ficticious magnetic field B& c oc J, the relative angular momentum. 
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Few-body problems have always attracted interest in the fields of atomic and nuclear 
physics. Recent work on laser-cooled ions in Paul traps [|IJ has heightened their importance 
in atomic physics. In condensed matter physics, such problems have been used indirectly as 
cluster calculations for understanding many-electron systems such as the two-dimensional 
(2D) electron gas in a magnetic field. A famous example is Laughlin's numerical calcula- 
tion for N = 3 electrons in a 2D parabolic potential used for investigating the Fractional 
Quantum Hall Effect ||. Few-body problems have recently taken on more direct relevance 
in semiconductor physics due to rapid advances in fabrication of quantum dots containing 
few electrons fl [|]- In lateral quantum dot structures, the electrons are typically 
free to move in only two spatial dimensions and the confining potential is approximately 
parabolic 0] [[|] . A complete description of this few-electron system is complicated since the 
confinement energy, the electron-electron repulsion and the cyclotron energy due to applied 
magnetic fields are typically comparable in magnitude. Numerical perturbative approaches 
employing a basis of non-interacting single-particle states become computationally intensive 
in the strongly-interacting (Wigner solid) regime. Analytic simplifications of the exact N- 
particle Hamiltonian or exact solutions of model iV-particle Hamiltonians can therefore be 
useful. 

Few-body Hamiltonians are rarely solvable analytically. Exceptions include N particles 
in ID with f3/r 2 interaction and N = 2 electrons in 2D with f3/r 2 interaction || and 
magnetic field. Here we show that the 2iV-dimensional problem of N particles (e.g. electrons) 
with (3/r 2 interaction in a 2D parabolic potential u (e.g. quantum dot) and magnetic field 
B reduces exactly to solving a (2N — 4)-dimensional problem which is independent of B and 
ujq. An exact set of relative mode excitations are obtained. The N = 3 particle problem 
reduces to that of a particle moving in a 2D non-linear potential of strength j3, subject to a 
ficticious magnetic field B^ c oc J, the total relative angular momentum. The ground state 
J (i.e. magic number) transitions for N = 3 are quantitatively consistent with numerical 
calculations for the Coulomb interaction |5[]. Analytic results are given in the Wigner solid 
regime. The present work implicitly includes mixing with all Landau levels. 
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The exact Schrodinger equation for N particles with repulsive interaction (3/r 2 moving in 
a 2D parabolic potential subject to a magnetic field B (symmetric gauge) along the z-axis, 
is given by (H space + H spin )^ = E^>\ 

^pace = + W" 2 0(B)\n\ 2 + y^) + E uT-^Tp W 

i=l " i<.j I 1 3 1 
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where uJq(B) = uol + c<j c is the cyclotron frequency, and -ff sp in = —g*I^BBJ2i s i,z- The 
momentum and position of the i'th particle are given by 2D vectors pi and respectively; 
U is the z-component of the angular momentum. The exact eigenstates are written in terms 
of products of spatial and spin eigenstates obtained from H space and H spin respectively; 
eigenstates of H spin are just products of the spinors of the individual particles. We employ 
standard Jacobi coordinates Xj (i — 0, 1, . . . , N — 1) where X = — J2j r j (center-of-mass), 
Xx = v /i(r 2 - X 2 = Jli^Y 2 ^ - r 3 ) etc. (see Fig. 1 for iV = 3) together with 
their conjugate momenta P«. The center-of-mass motion decouples, H space = HqmO^o) + 
iJ re i({Xj >0 }), hence E space = E CM + E re \. The exact eigenstates of -f^cM and energies E C m 
are well-known || . The non-trivial problem is to solve the relative motion equation H re \ifj = 
E re \ip. We transform the relative coordinates {Xj>o} to standard hyperspherical coordinates: 
Xi = riJJ^ 2 smaj + i)cosaie' l9i with r > and < on < | («i = 0). Because J remains a 
good quantum number, we introduce a Jacobi transformation of the relative motion angles 
in particular & = ^E^I 1 ^, Q = Q x -Q 2 etc. (see Fig. 1 for N = 3). We hence 
have (N — 1) variables, (N — 2) a- variables and a hyperradius r giving a total of 2(N — 1) 
variables as required for the relative motion. The exact eigenstates of if re i have the form 
i/j = e iJB ' R{r)G{VL) where denotes the (2N - 4) remaining {8, a} variables excluding 9'; 
R(r) and G(Q) are solutions of the hyperradial and (2iV — 4)-dimensional hyperangular 

equations respectively. The hyperradial equation is 

l_ 2N-3 d _ 7(7 + 2N - 4) _ r 2 2m*(ff rcl - hJf ) _ 
l dr 2 + r dr r 2 Z 4 + Al 2 ^ 

where Zq = ft(m*u;o(.B)) _1 ; the parameter 7 > and is related to the eigenvalue of the B 
and cuo- independent hyperangular equation (see below). Equation (2) can be solved exactly 
yielding 



E rel = kuj (B)(2n + ^ + N - I) + (3) 
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R{r) = { T VLl +N -\- I )e~^ . (4) 

In in 



where n is any positive integer or zero and 

; / -^n I 72 

Equation (3) provides an exact (and infinite) set of relative mode excitations 2huo(B)An 
for any iV regardless of particle statistics and/or spin states. These are "breathing" modes, 
as shown below for N = 3; numerical Coulomb results have shown similar modes to this set 



of /^-independent excitations |TIJ. All that remains is to solve the B and wn-independent 
hyperangular equation which resembles a (single-particle) Schrodinger-like equation in (2iV— 
4)-dimensional f2-space. The eigenvalue of the hyperangular equation 

e= ^ 7(7 + 2W _4)- J »-( s ^)1 (5, 

where Kiass is the potential energy of the classical, minimum-energy configuration (Wigner 
solid); Vdass oc (3~2Uq{B) and hence e (like 7) is independent of B and ujq. The exact relative 
energy for any iV 

E ml = hu (B)[2n + ([N - 2f + J 2 + (j^f + |) § + 1] + ^ • (6) 

^rei only depends on particle statistics through e. As % — > 0, E Te \ — > V c i ass and e — > 0. 
Physically, e accounts for the "zero point energy" in f2-space associated with the quantum- 
mechanical spread of G(Q) about the hyperangles Q corresponding to the classical, minimum 
energy configuration (Wigner solid); the actual spread in G(O) and hence e will depend 
on total wavefunction symmetry requirements (see below for N = 3). In general e > 0, 
e ~ (3^ where /i < 1 (the dominant /^-dependence of E Te [ lies in (Kiass) 2 ) an d e ~ J 5 where 
5 < 2. Equation (6) implies that for any N, the ground state J value will tend to become 
increasingly large and negative with increasing I?-field (u c > 0, e.g. electrons). We now 
demonstrate these statements explicitly for N = 3. 

For iV = 3 we change variables from a, 9 to x, y where x = ln(tana) and y = 9 — \. Since 
< a < f , hence —00 < x < 00 (N.B. —n<y< n). We define p x = |^ and p y = | 
The exact hyperrangular equation can be written in the form 



2 / I ftJcos(2tan 1 e x )\2 

where 

„, , *„ r (2 + cos(2tan~ 1 e x ) 3 . 2/fVL , 

War.uje) = m*/3 7 ; \ — r —, 9 -sm i (2tan V) 

v ' y ' ; PL (cosec(2tan"V) +cot(tan"V)) 2 -3sin 2 ?/ 4 v 1 

+-cos 2 (tan-V) + ^-cos 2 (2tan- 1 e :E )] . (8) 
2 m*p 

Equation (7) represents the single-body Hamiltonian for a ficticious particle of energy e and 
unit mass, moving in the xy-plane in a non-linear (i.e. e-dependent) potential V{x, y;e), 
subject to a ficticious, non-uniform magnetic field in the ^-direction 

5 fic = ^[l-cos(4(tan-V))] . (9) 

Bn c is independent of B and has a maximum of at x = for all y. For small x, 
Bfi c »s — ^ 2 )- As a; — > ±oo, _B fic — > 0. Note we have here chosen to highlight the 

Schrodinger-like form of Eq. (7); a simple rearrangement of Eq. (7) shows it to be hermitian 
with a weighting function sin 2 (2tan _1 e :c ) |]TT|. Our results are exact so far. 

Figure 2 shows the potential V(x, y;e) in the (x, y) plane. V(x, y;e) > everywhere. 
Minima occur at (0, 0) and (0, ±7r) where V(x, y; e) = (N.B. (0, tt) is equivalent to (0, — 7r)). 
Maxima occur at (ln\/3, ±|) in Fig. 2, where V(x, y; e) — ► oo. Since e > 0, these statements 
hold for any e. We now discuss the physical significance of these features. The classical 
configurations of minimum energy (Wigner solid) correspond to the particles lying on a ring 
in the form of an equilateral triangle with Ki ass = u (B)[6m* j3]^ . There are two distinct 
configurations with clockwise orderings (132) and (123) corresponding to (a, 9) = (|,±|). 
In (x,y) coordinates, these correspond to (0,0) and (0, tt) (equivalently, (0, — tt)). Hence the 
classical configurations coincide with the minima in V(x, y; e) in Fig. 2 and the maximum in 
Bfi c . The formation of the Wigner solid should therefore be favored by both large B^ c (i.e. 
large |J|) and deep V(x, y;e) minima (i.e. large ft, strong electron-electron interactions). 

Consider the limit of three electrons with very strong electron-electron interactions (i.e. 
(3 — > oo). Since the tunnel barrier height between the two V(x,y;e) minima ~ /3, the 



ficticious particle sits at one of these minima and the system is locked in one of the two 
classical configurations, e.g. (132) at (0,0). The tunnelling probability between the minima 
is zero. Tunnelling between the two minima implies a mixture of configuration (123) into 
(132) and hence interchange of the original electrons; in many-body language exchange 
effects arising from wavefunction antisymmetry are therefore negligible, e is small compared 
to m* (3 and Eq. (6) reduces to 



6m* (3 i T huJc 



E rel = hu (B)[2n + (l + J 2 + + • (10) 

The energy E ve \ > Kiass since it includes the hyperradial zero-point energy (N.B. h — > 
yields E vcl -> Kiass and B &c 0). 

Next consider large but finite f3. The ficticious particle now moves in the vicinity of the 
minimum (i.e. (x,y) ~ (0,0)). The electrons in the Wigner solid are effectively vibrating 
around their classical positions. Expanding the potential V(x, y;e) about (0,0) to third 
order, the exact Eq. (7) becomes 



[ 2 2 2 x 2~ v 



+ 2 7 + ~y x x 2 + -u 2 y y 2 ]G(x, y) = eG(x, y) (11) 



where uj 2 x = ( 3m A 13 + 2e) and uj 2 = 3 ™ 13 . This has the form of a single electron moving in 
an anisotropic parabolic potential, subject to a uniform magnetic field B^ c = Equation 
(11) is exactly solvable for e using a symmetric gauge |12[ (the energies are independent 
of the choice of gauge for Ba c ). A full discussion of the results for any e will be presented 
elsewhere. As an illustration, we consider small e hence u> x ~ u y . Equation (6) becomes 

E rel = hu (B)(2n + [1 + J 2 + + 2(2n' + \l\' + 1)(J 2 + ^S^)^ + 21' jf* + 1) 

■ (12) 

The ficticious particle has its own set of Fock-Darwin (and hence Landau) levels labelled 
by n' and a ficticious angular momentum I'. For large /3 and small n', I' and J, Eq. (12) 
yields an oscillator excitation spectrum with two characteristic frequencies y2hujQ(B) and 
2hu) (B) representing "shear" and "breathing" modes of the Wigner solid. 



For smaller (3 (i.e. weaker interactions) and/or larger e (i.e. excited states), the tun- 
neling probability between the V(x, y;e) minima in Fig. 2 becomes significant; the Wigner 
solid begins to melt and wavefunction antisymmetry (exchange) must be included. For three 
spin-polarized electrons, ip must be antisymmetric under particle interchange i <-> j. The 
hyperradial part R(r) is invariant; particle permutation operations in (i - !,^,^) become 
straightforward space-group operations in the (x,y) plane. For small (x,y), 1 «-> 2 is equiv- 
alent to (x, y) — > (x,y + it) with 6' — > 6' + | ; 1 «-> 3 is equivalent to (x, y) — > (x,y — it) with 
0' — >■ 0' + | ((x, represents (x, rotated by 4f); 2 «-> 3 is equivalent to (x, y) — > (x, y + 7r) 
with 6 1 ' — > 6' — | ((x,y) represents (x, y) rotated by — Single-valuedness of ^ implies 
e ±jj7r G(x, ± 27r) = G(x, y). Note we have implicitly satisfied Bloch's theorem in this anal- 
ysis, i.e. G(x,y ± 2n) = e ±l2nk G(x,y). The solutions G(x,y) of Eq. (7) with the lowest 
possible e (and hence lowest E Te \ at a given u c ) should be nodeless in the vicinity of (0, 0) (c.f. 
ground state in the parabolic potential with n' = = I' in Eq. (12)). However the above 

2 J 

symmetry requirements forbid such a nodeless solution unless e ln ~ = 1. Therefore the only 
symmetry- allowed solutions G(x,y) which are nodeless (i.e. smallest e and hence lowest E re i 
at a given uj c ) are those where J is a multiple of three. Evaluating the simplified expression 
for E Te \ in Eq. (12) (n' = 0= /'), the following ground state J transitions are obtained with 
increasing uo c for three spin-polarized electrons in a GaAs dot (Huq = 3.37meV as in Ref. 
@) H; -3 ->• -6 at B = 5.0T, -6 -> -9 at B = 8.7T, and -9 -»■ -12 at B = 12.2T 
(N.B. J = is not allowed by symmetry). The numerically obtained values from Ref. 
are B ~ 5.5T, 8.4T and 12. 4T using a Coulomb interaction. Our analytic results therefore 
agree well with the numerical calculations despite the different interaction form (see below). 
A feature of these analytic results is that they become more accurate in the Wigner solid 
regime (e.g. large f3 or \ J\) while the numerical calculations become more computationally 
demanding. 

For general N, the hyperrangular equation (c.f. Eq. (7)) becomes (2iV — 4)-dimensional. 
However in the Wigner solid regime (large (3 or |J|) the classical minimum energy config- 
urations will still be important in determining e and hence E Te \, just as for N = 3. The 

7 



classical minimum energy configurations (with 1/r interaction) for iV < 6 all consist of iV 
particles on a ring; for N = Q additional minima occur ||14|| . Intriguingly it is at iV = 6 
that the magic number J sequence of A J = N is broken jHJ. The present formalism which 
emphasises classical configurations may shed light on a possible link here. 

Finally we note that the f3/r 2 interaction (f3 > 0) is not unrealistic in quantum dots 
due to the presence of image charges; in particular it resembles the dipole-like form used 
successfully in Ref. 4. Furthermore, recent quantitative comparisons fl2|] ]10| have 
shown that the 1/r 2 and 1/r repulsive interactions yield energy spectra with very similar 
features (e.g. ground state J transitions, the relative excitation 2hu>o(B) for N = 2 [|Kj); 
the above results for N = 3 are consistent with this finding. Significant differences will only 
arise for the case of attractive forces f3 < (e.g. between electrons and holes) because of 
the increased importance of the r — > dynamics for that case. 

This work was supported by COLCIENCIAS Project No. 1204-05-264-94. 
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Figure Captions 



Figure 1: The N = 3 system. Reading clockwise, classical configuration for three repulsive 
particles (132) corresponds to (a, 9) = (f,§) (i.e. (x, y) = (0,0)); (123) corresponds to 
(a, 9) = (f,-f) (i.e. (x,y) = (0, tt) or, equivalently, (0,-tt)). 

Figure 2: Contour plot of ficticious potential V(x, y;e) in the (x,y) plane for N = 3. 
Relevant corresponding configurations are shown. Minima in V(x,y;e) occur at (0,0) and 
(0, ±7r) (i.e. at classical configurations). Maxima occur at (ln-\/3, ±§ ), where V(x, y;e)^oo 
(i.e. particles 2 and 3 or 1 and 3 coincident). V(x, y;e) is positive and finite everywhere 
else. The same qualitative features appear for all e {-^rs = 5 for illustration). 
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